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I. INTRODUCTION 



In the last few years the magnetic catalysis (MC) of chiral symmetry breaking [1]- [3] has been the focus of attention 
of many works on non-perturbative effects of magnetic fields [1]- [21]. The phenomenon consists on the dynamical 
generation of a fermion condensate (and consequently of a fermion mass) when the fermion interactions occur in the 
• presence of an external constant magnetic field. A most significant feature of the MC is that it requires no critical 
C*") ' value of the fermion's coupling for the condensate to be generated. That is, the symmetry breaking takes place at the 
ON ! weakest attractive interaction. Physically it is due to the fact that the magnetic field forces the low energy fermions 
to reside basically in their lowest Landau level (LLL), while the higher energy fermions actually decouple [8]. This, in 
turn, yields a dimensional reduction of the infrared fermion dynamics. The dimensional reduction is reflected in an 
i effective strengthening of the fermion interactions leading to dynamical symmetry breaking through the generation 
i-Q ' of a fermion condensate. 

A particularly important question to understand in this context is how the MC is affected by the introduction of 
' fermion-scalar interactions. Fermion-scalar interactions are an essential element of the unified theories of fundamental 
forces. As is well known, they are expected to be responsible for the fermion mass appearing due to the spontaneous 
i-C ■ symmetry breaking of the electroweak symmetry. Fermion-scalar interactions are also relevant in condensed matter 
| physics where the complexity of strongly-correlated many-body systems some times calls for a description in terms of 
more simple, phenomenological theories that contain interacting scalar in addition to fermions (see e.g. [22]). 

In Refs. [12]- [13] two of us studied the realization of magnetic catalysis in a (3+l)-dimcnsional Higgs- Yukawa (HY) 
model, showing that the magnetic-field-induced fermion mass is enhanced by fermion-scalar interactions. As we will 
show below, this enhancement is also found within a more accurate approximation for a wide range of couplings. This 
result might find applications in early universe transitions, as well as in condensed matter physics. 

In [12]- [14] some applications of the MC to the early universe were briefly considered. They were motivated by 
many astrophysical observations of galactic and intergalactic magnetic fields indicating the existence of seed fields 
that originated from large primordial magnetic fields (for a recent review on cosmic magnetic fields see [23]). If the 
primordial magnetic fields in the early universe were large compared to the values close to the phase transition point 
of the fermion masses generated through the usual mechanism of spontaneous symmetry breaking, the fermion would 
seem approximately massless. Under these circumstances, it is important to investigate if the primordial magnetic 
fields could contribute to the masses of the fermions through MC and hence influence the phenomenology of the early 
universe [13]. 

On the other hand, to discuss applications of MC in the context of a HY theory to condensed matter, we need, 
besides interactions modelled by fermion-scalar terms, a physical system that, despite being non-relativistic, can be 
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described under certain conditions by a " relativistic" Hamiltonian. We will see below that these conditions are indeed 
present in the physics of high-T c superconductors. 

High-T c superconductors, which are characterized by the existence of nodal points where the order parameter (gap 
function) vanishes, provide a practical realization of a "relativistic" system in condensed matter physics. This is so 
because the low-energy spectrum of the nodal quasiparticlcs is linear, hence the quasiparticle excitations are described 
by an anisotropic Dirac Hamiltonian [24]. In Rcf. [22] a quantum-critical phase transition to a new superconducting 
state, characterized by the appearance of a secondary pairing at some doping level, was proposed to explain recent 
measurements [25] of an anomalously large inelastic scattering of quasiparticles near the gap nodes of a superconductor. 
The observed secondary pairing transition made the nodal quasiparticlcs fully gapped. Based on the symmetries of the 
superconductor, the authors of Ref. [22] made a classification of a set of fermion-scalar interactions that in principle 
could be in agreement with the experimental observations, and then performed a perturbative renormalization-group 
analysis of each model to determine the possible existence of a quantum-critical point. In Ref. [26], expanding on 
the ideas of [22], the existence of a quantum critical point was established directly in a (2+l)-dimcnsional HY theory 
beyond a non-perturbative approach, which allowed to make quantitative predictions for the corresponding quantum- 
critical behavior. The gap generation (fermion mass) was associated in [26] with the breaking of a discrete chiral 
symmetry. 

We would like to underline that the breaking of the chiral symmetry in [26] was found to occur when the Yukawa 
coupling (assumed to be related to the doping level) reached a critical value, that is, the symmetry breaking was 
not associated to the phenomenon of MC, as no external magnetic field was introduced in the analysis. However, 
as recently observed [27] by measuring the splitting of the conductance peak that characterizes the nodes of high-T c 
superconductors, the development of a secondary quasiparticle gap may be triggered not only by the doping level, but 
also by an applied magnetic field. Could the secondary gap triggered by the magnetic field be the consequence of MC 
occurring within the superconductor? We believe that the results we are going to derive below strongly indicate that 
the answer is yes, if, as argued in [22] and [26], the HY theory is the model describing the appearing of the secondary 
gap. Nevertheless, to match the experimental observations we would need to particularize the analysis done in the 
present paper to the (2+l)-dimensional case, and adjust the physical values of the couplings to those characteristic 
of the superconductor. 

As already mentioned, in Rcf. [13] the phenomenon of MC in a (3+l)-dimcnsional Abclian gauge theory with HY 
interactions was studied. In that work it was shown that the non-perturbative solution of the minimum equations for 
the composite-operator effective action leads not only to a magnetically catalyzed fermion dynamical mass, but also 
to a nonzero scalar vev ip c and consequently, to a nonzero scalar mass. In other words, thanks to the magnetic field, 
a scalar-field minimum solution is generated by non-perturbative radiative corrections. 

We should underline though that the fermion and scalar masses of rcf. [13] were obtained within a simplified 
approximation known in the literature as the constant mass approximation (CMA) . In general, to find the dynamical 
mass — which is nothing but the part of the fermion self-energy proportional to the identity matrix — one has to solve 
a non-perturbative gap equation (i.e. the Schwinger-Dyson equation for the full fermion propagator). This means to 
solve a non-linear, implicit integral equation for the fermion self-energy, which is a momentum-dependent function. 
Most authors approach such a mathematically complicated problem with the help of the rough CMA approach. It 
consists on neglecting the momentum dependence of the self-energy in the gap equation. This is done by substituting 
the self-energy function in the gap equation by its value at zero momentum, that is, by the infrared mass. There is 
no general principle that guarantees the validity of this approximation for the whole range of physical couplings. 

For theories with several couplings, due to the richness of the parameter space, the reliability of the CMA is 
questionable and should be investigated in detail. In the case of the HY model, aside from the multiple-coupling 
problem, one has to deal with a system of non-linear, coupled integral equations, one for the fermion dynamical mass 
and other for the scalar vev [13]. We cannot disregard in this situation the possibility of regions of these parameters 
where the CMA is reliable and regions where it is not. In this case one has to turn to a more accurate approximation 
on which the momentum dependence of the self-energy is taken into account when solving the gap equation. This 
more accurate approximation is known as the beyond-constant-mass approximation (BCMA). 

For theories like QED containing only one coupling constant, the CMA is known to be appropriate, since going 
beyond it does not produce qualitatively different results. This has been explicitly shown for (3+l)-dimcnsional [2] 
and (2+l)-dimensional QED [28], and independently corroborated by our calculations below. In Rcf. [9] the BCMA 
mass solution of (3+l)-dimcnsional QED was found to agree with the CMA mass obtained from the improved-ladder 1 
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gap equation. This result was later corroborated by numerical calculations in [29] . 

Considering that different physical applications of the MC in models with fcrmion-scalar interactions would require 
different values of the couplings constants, and in particular, given the relevance that the Abelian gauge Higgs- 
Yukawa theory may have for condensed matter and other field theory applications, it is important to perform a 
BCMA investigation of this model in all possible regions of the parameter space and find out whether it significantly 
differs or not from the CMA results. A main goal of the present paper is to carry out such a study. 

By going beyond the CMA, we will determine the region of Yukawa and scalar self-interaction couplings where the 
CMA is valid, and will obtain the numerical BCMA solutions for the fermions and scalar dynamical masses in the 
complete physically meaningful parameter region. As we will see below, the CMA results for the Abelian gauge Higgs- 
Yukawa theory are mostly reliable in the available parameter space. An important finding is that the (BCMA-found) 
mass values are many orders of magnitude larger than those obtained in the absence of fermion-scalar interactions, 
corroborating, within this more accurate approximation, the enhancement of the dynamical mass by the Yukawa term. 

The paper is organized as follows: In Section II we derive the non-linear integral equations for the fermion self 
energy (gap equation) and the scalar vev in a gauge Higgs- Yukawa theory. The integral gap equation is then converted 
into a second order differential equation with boundary conditions. In Section III, this boundary-value problem is 
analytically solved leading to the self energy as a function of the momentum and the infrared fermion mass. Using 
the self-energy solution and the equation for the scalar minimum, we arrive at two coupled transcendental equations 
depending on the infrared dynamical mass and the scalar vev. These equations are numerically solved and the results 
are used to determine the region of reliability of the CMA and to compare the mass values obtained in the CMA and 
in the BCMA approaches. We end Section III discussing the solution of the gap equation at zero Yukawa coupling, 
and showing that it leads to the same result found in Ref. [2] for (3+l)-dimensional QED. In Section IV, we state our 
concluding remarks and reconsider the question of the relevance of the magnetic catalysis in the electroweak phase 
transition using the BCMA results. 



II. INTEGRAL EQUATIONS 

Let us consider the following Lagrangian density 

l = -\f» v f„ v + + e^yvA, - \d^d^ - ±y - - (i) 

that describes a gauge Higgs- Yukawa model with a fermion field coupled to scalar and electromagnetic fields. The 
scalar field is electrically neutral, but self-interacting. 
The Lagrangian density (1) has U(l) gauge symmetry, 

An -> + ^df,a(x) 
*l> -» e ia{x) ip, (2) 

fermion number global symmetry 

4> - e l >, (3) 

and discrete chiral symmetry 

i\) — > 75^; — > — V>7 5 , <f — P- (4) 

Notice that the quadratic scalar term has the correct sign of a mass term, thus no vacuum expectation value of 
the scalar field exists at tree level. In the course of our calculations we will take /x — > to search for a dynamically 
induced mass. The discrete symmetry (4) forbids a mass for the fermions to all orders in perturbation theory. 
Nevertheless, this symmetry could be dynamically broken through non-perturbative generation of a composite field 
(fermion-antifermion condensate). Such a fermion condensate would lead to a dynamical fermion mass and to a 
non-zero vacuum expectation value of the scalar field [13], which in turn would contribute to the scalar mass. 



the mass solution (either CMA or BCMA), assuming that the gap equation is found using a ladder or an improved ladder 
approximation. 



3 



It is known that in the case of the non-gauge (3+l)-dimensional Higgs- Yukawa theory, no value exists for a running 
\ y at which a chiral symmetry breaking fermion condensate can be generated 1 . As shown in [13], the situation 
drastically changes when a magnetic field is introduced. In this case a non-trivial solution exists at the weakest value 
of \ y and one can show that a fermion condensate (together with a dynamical fermion mass and a scalar vev) is 
magnetically catalyzed. 

However, as already mentioned, the solutions in Rcf. [13] were found within the CMA, and therefore it is important 
to investigate their reliability beyond that approximation. Our task hereafter will be to extend the results of Ref. [13] 
beyond the CMA to find the dynamical mass and the scalar vev for all physically meaningful values of X y and A. For 
the sake of understanding, we will repeat the outline of the derivations done in Ref. [13] that lead to the coupled set 
of integral equations (gap and scalar vev equations) that will be the starting point of our new calculations. 

Let us consider the Lagrangian density (1) in the presence of an external constant magnetic field B (without loss 
of generality we assume that the magnetic field is directed along the third coordinate axis and that sgn (eB) > 0), 
which can be introduced by adding the external potential = (0, 0, eBxi, 0) as a shift to the oscillatory gauge field 
in Eq. (1). To find the vacuum solutions of this theory we need to solve the extremum equations of the effective 
action T for composite operators [30], [31] 



8G 



= 0, (5) 
= (6) 



In the above G(x,x) — a(x) — (0 | vp(x)ip(x) | 0) is a composite fermion-antifermion field, and ip c represents the vev 
of the scalar field. The subindex B indicates that the effective action is considered in the background of the external 
magnetic field. 

Equations (5) and (6) are, respectively, the Schwinger-Dyson (SD) equation for the fermion self-energy operator E 
(gap equation) and the minimum equation for the vev of the scalar field. As we are interested in the possibility of a 
scalar mass induced — through the interactions with the fermions — by a dynamically generated fermion condensate, 
we will set, as stated above, the bare scalar mass [i to zero. Notice that, if the minimum solutions of Eqs. (5) and (6) 
are non trivial, the discrete chiral symmetry (4) is dynamically broken and both fermions and scalars acquire mass. 
The loop expansion of the effective action T for composite operators [30], [31] can be expressed as 

r B (G, ifc) = S (<p c ) - iTr In + i-Tr \iiD~ 1 + i-Tr In A" 1 

-iTr[G- 1 (^c)G\+T 2 (G^c)+C (7) 

Here C is a constant and S (ip c ) is the classical action evaluated in the scalar vev ip c . Non-bar notation indicates free 
propagators, as it is the case for the gauge 

and the scalar 

f d 4 a e iq-(x-x') 

propagators. Here £ is the gauge fixing parameter and M 2 = ^ip 2 denotes the scalar square mass. A dependence on 
full boson propagators is not included since we do not expect the gauge field to acquire nonzero expectation values 
for its composite operator. On the other hand, we are going to explore the possibility of a non-zero vev of the scalar 
field, hence, a composite-operator solution for the scalar would be a correction of higher order that can be neglected. 

The bar on the fermion propagator G (x, y) means that it is taken full. The full fermion propagator in the presence 
of a constant magnetic field B can be written as [5], [11], [32]- [33], 



1 The incorporation of gauge field terms in the Higgs- Yukawa model may lead to chiral symmetry breaking at some critical a, 
just as it occurs in (3+l)-QED [1]- [2]. 
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with S(p) being the fermion self energy, p = (po, 0, — ^/2gBk,ps), and fc denoting the Landau level number. Similarly, 
the free fermion inverse propagator in the presence of B is given by 

G- 1 (x, y, y> c ) = W dP0 ^ P2d / 3 E p (x) (j.p + \ y <p c ) E p (y) . (11) 
k J ( 2?r ) 

Note that X y <p c enters as a contribution to the fermion mass due to the shift ip — > tp + ip c in the scalar field done in 
the classical action to account for a possible non-zero scalar vev. The value of <p c will be determined self-consistently 
through Eq. (6). 

In the above equations, Ritus' E p functions [32]- [33] were introduced. They form an orthonormal and complete 
set of matrix functions and provide an alternative method to the Schwinger's approach to problems of QFT on 
electromagnetic backgrounds 2 . Ritus' approach was originaly developed for spin-1/2 charged particles [32]- [33], and 
it has been recently extended to the spin-1 charged particle case [34]. 

The function T 2 (G, <p c ) in (7) represents the sum of two- and higher- loop two-particle irreducible vacuum diagrams 
with respect to fermion lines. For weakly coupling theories, like the case of Lagrangian (1), one can use the Hartree- 
Fock approximation, which means to retain only the contributions to T 2 that are lowest-order in coupling constants 
(i.e. two- loop graphs only), so that it becomes 

r 2 (G,^ c ) =y J d 4 xd 4 ytr[G(x,y)^G(y,x)^D^(x,y)} 

J d 4 xd 4 ytr (^G{x,x))D^{x-y)tr [fU{y,y)) 

Xy - f d 4 xd 4 ytr[G(x,y)G(y,x)A(x,y)] 



2 

Ay 



J d 4 xd A ytr (G (x, x)) A(x - y)tr (G (y, y)) (12) 



As discussed above, the infrared dynamics (p « \j2eB) of a system of interacting fermions in the presence of a 
magnetic field is mainly governed by the contribution of the LLL [1]- [2]. To obtain an explicit form for Eqs. (5)- (6), 
we use the propagators (8)-(10) in Eqs. (7) and (12), and take into account that in the background magnetic field 
the self-energy structure entering in the full fermion propagator (10) should be written as [11] 

E(p) = Z„ (p) 7 .p B + Z ± (p)j.p ± + S(p). (13) 

Here we are using the notation p = (po,P3) and p± = (pi,p 2 ) for the momentum components. The wave function 
renormalization coefficients Z„ ,Z±axe scalar functions of the momentum. Using this structure for £ in the full fermion 
propagator, evaluating at the LLL (k=0), and using the solution of the wave function renormalization, = 0, found 
in Ref. [12], we have that the gap equation (5) and the scalar minimum equation (6) of our theory take the form 

2 ]~d\d% S((«-p)J) e -?i °fd%diq i S((g-p) 2 ) e -5i 



^(p) = 2^ / ^ - ^ 2 , ~-± ~ + y 



and 



(2nf {q-p) 2 + ^((q-p) 2 ) <Z 2 y J (2nf (q- p)J + £ 2 ((? - §2 + M 2 

-p-J^ d% ^_ + A^ c (14) 



0t = —JL- / d 2 q„ ^ (15) 

2AW ^ + £ 2 (^) 



2 For details on the use of Ritus' method in the theory given by Eq. (1) see Refs. [12]- [13]. 
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respectively. Dimensionless field- normalized quantities are denoted by Q = S = . Notice that if we set X v — in the 

above equations, Eq. (14) reduces to the same gap equation found in [2] for (3+l)-dimensional QED, since, in the 
absence of a Yukawa term, the theory (1) becomes equivalent to a QED theory on which an extra, but disconnected, 
real scalar field has been added. 

Changing q n to polar coordinates (k, 8) in the above integrals and integrating in the angle, we find 

£(p) = SL [ d fc2 SCfc 2 ) 2 ~ 2 /■ ~ 2 ^(fc 2 ) ( g-2 fe 2) 



OO 



4tt 2 m 2 J 



dk z ^ v ^ ^ + A y £ c , (16) 

fc 2 + £ 2 (fc 2 ) y 



ft 



oc 

2Xtt 2 J 



dfc 2 



S(fc 2 ) 
fc 2 + £ 2 (fc 2 ) 



(17) 



The functions >rt(p 2 ,2:) are defined by 



^t(p 2 ,z) 



<iz- 



+ p 2 + x + t) 2 - ixp 2 



(18) 



To make the calculation more manageable, it is convenient to divide the momentum integration in Eq. (16) in the 
two regions separated by the dimensionless momentum square p 2 . Expanding the kernels *c t (p 2 ,x) appropriately on 
each region, we find 



( -2 

V 



m - s < 



dk 2 



E(fc 2 



k 2 + £ 2 (/c 2 ) 



p 2 + z 



+ / & 2 



£(fc 2 



+ 



16tt 2 



v 

1 



dk 2 



£(fc 2 ) 



/ 



dz- 



k 2 + Y, 2 (k 2 )J p^ + z + M 2 

p 



k 2 + E 2 (k 2 ) 

OO 

+ J dk 2 



dz- 



k 2 + . 



S(fc 2 



dz— t^— 

k 2 + Y, 2 (k 2 )J k 2 + z + M 2 



2 ^ - 



(19) 



Notice that we used Eq. (17) to combine the last two terms of Eq. (16) into the last term of Eq. (19). The analytical 
solutions of Eqs. (17) and Eq. (19) can be explored by converting first the non-linear integral equation (19) to a 
second order non-linear differential equation. First, however, we must take into account that the consistency of the 
LLL approximation requires to use a momentum cutoff of order ^/2eB in the momentum integrations, and hence the 
infinity limit in all the integrals in k 2 should be changed to 1. 

One can easily see, by taking derivatives of Eq. (19) with respect to x = p 2 and combining them conveniently, that 
the integral equation (19) is equivalent to the following second order differential equation 



t"{x)-%^-t'{x)-g'{x) ^ {X) 







g'(x) v ' ^ ' x + s 2 (x) 

If we now differentiate (19) and evaluate the result at x = 0, we obtain the following boundary condition 

E'(x) 



(x) 



0. 



x=0 



where 



9(x) = ^- g {x) + ^g(x + M 2 ), 



(20) 



(21) 



(22) 
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9(V) = J dz^- y (23) 
o 

Similarly, taking the derivative of (19), multiplying it by JtJ^ and evaluating at x = 1, we obtain the second 
independent boundary condition 

2 

= ^yVc- (24) 

x=l J 

In doing so, we have traded a non-linear integral equation for a non-linear boundary value problem. Finding the 
solutions to the coupled set of Eqs. (17) and (20), with boundary conditions (21) and (24) will be the aim of the next 
section. 
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III. FERMION AND SCALAR MASSES IN THE BEYOND-CONSTANT-MASS APPROXIMATION 

A. Beyond-Constant-Mass Analytical Solutions 

An analytical expression for the solution E(x) of (20)-(21) can be found considering a linearized version of the 
equations (20) and (17), on which the fermion self energy in the denominators is replaced by its zero momentum value 
£ (0) = m. The consistency of such linearization is justified if the self energy is a rapidly decreasing function of the 
momentum. We will corroborate at the end of the derivations that follow that this is indeed the case. Then, the gap 
equation (20) can be written as 



while the equation for the scalar minimum takes the form 

i ^ 

$l = ^\l dx^l^ (26) 
2Xn 2 J x + m 2 

o 

From a physical point of view, we expect that the masses for both fermion and scalar fields will be much smaller 
than the magnetic field that induces them through the formation of a fermion condensate. Therefore, it is reasonable 
to assume that m 2 <C 1 and M 2 < 1. At the end of our calculations we must check in the obtained results the 
consistency of this assumption. 

Taking into account the asymptotic behaviors of the function g (x) in the regions: 



(1) x « M 2 « 1 



(2) M 2 «x<l 



g'(x) x y w 2ttx v ' 



1 

g' (x) x 




H^-k + T^ - (28) 



one ends up with a different boundary value problem at each region. The two boundary value problems are defined 
by the following equations: 



(1) For x « M 2 « 1 



E" (x) + is' (*) + f- = 0, (29) 

X Z7T X (X + 771 ) 
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(2) For M 2 « x < 1 



xt' (x) x=Q = 



(30) 



a; \ Z7r I67H / x(x + m^) 



= 0, 



(31) 



S(l) + eS'(l) = ?A^ c (32) 

where e = =^ = 1.477 and a = is the fine-structure constant. 

For most physically interesting applications of the gauge Yukawa theory, the Yukawa coupling X y is < 10 _1 . For 

those X'yS, the parameters v = and V = \J jr^ + practically coincide (for X y = 10~ 2 they already have three 
significant common figures). Thus, we can take V ~ v, in Eq. (31), reducing the problem to a single second order 
differential equation. The new problem is then defined by Eq. (29) and the two boundary conditions (30) and (32). 
The solution to this boundary value problem can be written as the following combination of hypergeometric functions 
(for properties and formulas of the hypergeometric functions see [35]) 

S (x) = A\F (iv, -iv; 1; -iL) + A 2 (l + ^) F (l - iv, 1 + iv; 2; 1 + ^ ) (33) 

Taking into account the boundary condition (30) and the formula 

dF (a,b:c; z) ab . _ , „ „ . ,„,. 
V ' ' ' = — F (a + l,fe+l;c+l;z) (34) 

we obtain A 2 = 0. As fh = E (0) , it is clear that A\ = fh. Therefore the self-energy solution becomes 

X (x) = fhF ( iv, — iv; 1; — ) (35) 

The second boundary condition (32) gives rise to 

1 \ v 2 ( 1 \ 2 

^2 - e — F 1 + *"» 1 - 2 ; -— = 0- 

TO Z / TO \ TO / 3 



mF it/, -it/; 1; -— ) - e—F I 1 + iv, 1 - iz/; 2; -— ) = -X v tp c , (36) 



which establishes a relation between the fermion dynamical mass fh, and the scalar vev {p c . This is an implicit, quite 
non-trivial equation for fh : besides the dependence on fh in the hypergeometric functions, the scalar vev (p c depends 
on fh through Eq. (26). 

To find the solution to the system formed by (26) and (36), we first note that Eq. (29) can be rewritten in the form 

±( x Y!{x))=-S*gL, (37) 
ax V / x -\- 7ji 



hence 

1 



dx , =-4 g, ( 1 )- (38) 

x + m z {x) v 1 



Using (33) and (38) in Eq. (26), and the values of A\ and A2 just found, we obtain 

^ = |^1f(i + ;m-^;2;-4) (39) 

From the asymptotic behavior of the hypergeometric function for large values of its argument [35] 

F (a b- c z) ~ r(c)F(6 - a) (-*)-' + r ( C )r(a-&) (_.r b (40) 
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we can show that 



F I iv, -iv; 1; 



T {-2iv) 



T {—iv)T (1 — iv) \m 



+ 



T (2iv) 



T(iv)T(l + iv) \m 2 



r(i + 2w) / 1 

r 2 (1 + iv) \fh 2 



T(l + 2iv) 



where 



and 



8 (v) = arg 



T 2 (l + iv 
T(l + 2iv) 



r 2 (i + w) 

r(l + 2ii/ 



^ 3 + O (i/ 5 ) 



r 2 (i + ^) 



so the function can be approximated by 

F ( iv, —iv; 1; 

Similarly, one can see that 



~ Re 



iu In 



= cos [ v In 



F 1 + 1- it/; 2;-— ~ sin [v In — 



1 



772 



(41) 



(42) 



(43) 



(44) 



(45) 



Substituting with Eqs. (44) and (45) into (36), (39), we obtain a much more simplified, though still transcendental, 
pair of coupled equations for the fermion infrared mass fh and the scalar vev (or equivalently, for the fermion infrared 
mass and the scalar mass M), 



t 4 A 4 

e~~[cos(t) - i/esin(t)] 3 - --^-sin(t) = 

9 ir z Av 



(46) 



^ = 2 = e ~ - -^(t)] 2 , (47) 

where the parameter t = i/ln(^j). 

Eqs. (46)-(47) represents the BCMA implicit solution for the fermion and scalar masses catalyzed by the magnetic 
field. This is as far as we can stretch our analytical calculations for m 2 and M 2 without introducing any additional 
approximation. In the following subsections we will perform a numerical analysis of these solutions. 



B. Numerical Solutions in the BCMA 



Since Eqs. (46)-(47) are highly transcendental, to obtain the explicit dependence with the couplings of the BCMA 
fermion and scalar masses, we have to resort to numerical methods. 

Figs. 1 and 2 display logarithmic plots of the numerical solutions of Eqs. (46), (47), versus couplings X y and A. 

From them, one can easily see that the two masses widely agree with the initial assumptions m 2 <C 1, M 2 <C 1. Only 
in the region of very large A y (very large n) and very small A (very small k) the fermion mass becomes of order one, 
hence, to be consistent, we should disregard the results in this corner. In any place out of this limited section of the 
parameter space, the results are reliable for both masses. 
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40 

FIG. 1. Logarithmic plot of the BCM solution field-normalized fermion mass square fh 2 versus couplings \ y = 10~ 6+ ~s and 
A = 10 _9+ 8 for a = — j=. The numbers in the x-y axes indicate the values of n and k, respectively. 




FIG. 2. Logarithmic plot of the dimensionless scalar square mass M 2 versus Yukawa and scalar self-interaction couplings. 
The intervals of the coupling constants used to generate this graph are the same as in Fig.l. 

Notice that the fermion mass grows with X y at any given value of A. This in turn implies an enhancement of the 
fermion mass as compared to its value within QED. While in QED the largest mass was no more than ~ 10~ 10 yf2eB 
[2], here the mass surpass this value in the majority of the parameter space in at least 5 orders of magnitude. 

It is because of such a significant enhancement of the dynamically generated mass in the presence of scalars, that 
the magnetic catalysis could play an important role in realistic applications of the HY model. The region of large 
\ y , large A, where the results are quite reliable, is the most interesting for applications to the electroweak theory, 
since the values of the coupling constants in that section include the value of the scalar self-couplings consistent with 
current experimental limits for the Higgs mass, as well as the Yukawa coupling of the top quark. 

To finish this subsection, let us consider the behavior of the self energy with the momentum. In Fig. 3 we have 
plotted the self energy solution (35) as a function of the momentum for fixed values of the couplings. As can be 
seen, decreases very quickly with the momentum. This behavior is in good agreement with the linearization used 
in Eq.(25). It also justifies the ultraviolet cutoff at y/2eB that was imposed on the integrals appearing in the gap 
equation (19) and the scalar minimum (17), since, as seen here, the main contribution to the integrals comes from 
the deep infrared region. 
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0.2 0.4 0.6 0.8 1 

FIG. 3. Self energy versus momentum squared for X y — 1CP 1 and A = 10~ 2 

C. Comparison Between the BCMA and the CMA Solutions 



To find the region of reliability of the CMA we have to determine the values of the couplings for which the two 
approximations give rise to the same mass solutions. With this aim we compare the BCMA equations (46)-(47) with 
the corresponding CMA equations 



TO = — A 1 ,lfi r H TO m —-rr 

3 y ^ 2 to 



(48) 



! -P'c 



3A y ^. 1 
-mm — , 

TTl 



2tt 2 A 



(49) 



that were previously found 3 in Ref. [13]. 

Eqs. (48)-(49) look very different from their BCMA counterpart: Eqs. (46)-(47). There is no reason to anticipate 
that the solutions of both sets of equations will coincide in all the parameter space. Combining Eqs. (48)-(49) we 
obtain 



9Att 2 
4A 4 



1-2, 2 1 

—vln—^ 



(50) 



From (50) we see that since to 2 has to be positive, the consistency of the CMA solution requires ^v 2 ln 2 ^j < 1> 
which is equivalent to have t < 1.4. Below, we will numerically check that this condition is indeed always satisfied. 

To compare the BCMA and CMA solutions we will explore whether there is a condition under which the CMA and 
BCMA equations reduce to an identical set. To this end, let us assume that VlnJ^ — vln^o <C 1. This restriction 

^ 5 Tin n\ A 

allows us to write Eqs. (46)-(47) as 



9^ 2 A 



(51) 



T?2 9A _ 2 
= 8A| m 



(52) 



B Here we have corrected some missprints appearing in Ref. [13]. 
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respectively. They are exactly the same equations found from (49) and (50), after using t<l. Thus, in this limiting 
case, the BCMA reduces to the CMA, thereby t <C 1 defines a condition of reliability of the CMA. 

The explicit region of parameter space where the CMA is reliable can be determined from a numerical plot of the 
ratio between the CMA and BCMA mass square solutions. To be sure that we are working with consistent masses, we 
will restrict the couplings to a strip in the (X yi X) plane, leaving out the corner of Fig. 1, where, as discussed above, 
the consistency of the approximation breaks down. 



log CM_BCM fermion mass square ratio 




log CM_BCM scalar mass square ratio 




FIG. 5. Log( + 2 CM ) in the region of couplings 10" 8 < A < 10"\ 10" 6 < X y < 10" 2 . 

M BCM 

Figs. 4 and 5 show logarithmic plots of the ratio of CMA over BCMA mass square results for fermion and scalar 
masses respectively, taken in the region of couplings 1CP 8 < A < 10 _1 , 10~ 6 < Ay < 10~ 2 . Both figures display 
similar behavior of the ratios, characterized by a discernible region of the parameter space, approximately given by 
10~ 4 < A < 1CP 1 and 10~ 6 < X y < 10 -5 , where a disagreement between BCMA and CMA results is apparent. 
However, even in this segment, the BCMA and CMA mass squares differ at most in one order of magnitude. Out of 
this limited region we find very good agreement between BCMA and CMA results, particularly at large X y , indicating 
that this is the most reliable region of the CMA solution within this model. 



12 



BCM - t - value 




CM - t - value 




The above observations are corroborated by the plots of the BCM- and CM- t's, as shown in Figs. 6 and 7 
respectively. Both surfaces have similar t-values at equal set of couplings, even when t 1 is not satisfied, indicating 
that, after all, and as already seen in Figs. 4 and 5, the two approximations give rise to very near mass values. Notice 
that the larger the Ay's, the smaller the i's in both approximation, leading to a better agreement between BCMA and 
CMA results, as expected from our previous analytical considerations. Therefore, although the numerical calculations 
show that the CMA results are widely reliable, it is in this extreme section of the parameter space where the two 
approximations totally coincide. From Fig. 7 it is evident that the CM-t never goes over the limiting value of 1.4, so, 
even in the region of larger discrepancy between the CM and BCM results (large A, relatively low Aj,), the CM mass 
solution remains real, as it should. The curves reflect the fact that the CM-approximation tends to overestimate the 
mass, because it substitutes in the integrals the self-energy function, which rapidly decreases with momentum, by a 
constant. 



D. BCMA in the X y = Limit (QED case) 

We shall discuss now the limiting case X y = which reduces to (3+l)-dimensional QED with a decoupled self- 
interacting scalar field. Let us find the solutions for the masses in this case. It is clear from Eq. (17) that no scalar 
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vev, and hence no scalar mass, is generated in this case. The fermion dynamical mass solution can be found from 
(46) evaluated at X y — 0. It leads to 

tan(t) = — . (53) 

In terms of ra 2 , it can be rewritten as follows 

m 2 = e -7T arctan (777) (54) 

Taking into account that ve <C 1 and using the asymptotic behavior arctan(x) ~ ^ — i, we obtain 

m 2 = e e e~^^- (55) 

This result coincides with the BCMA results found for QED within the ladder approximation (see Ref. [2] for details). 
As known, it is qualitatively very close to its CMA counterpart to 2 ~ e _7r ^° [2]. Thus, we are corroborating here the 
conclusion of the authors of Ref. [2], namely, the reliability of the CMA approach in (3+l)-QED 4 . 

It is worth to notice that the dynamical mass behavior is basically affected by the infrared conditions of the self 
energy, but it is practically indifferent to the ultraviolet boundary condition used in Ref. [2]. This explains why, 
despite using a momentum cutoff at ^/2eB and imposing the second boundary condition at x = 1, we still get in the 
\ y = case the same result as in [2] , where the momentum was allowed to run up to infinity. 



IV. CONCLUDING REMARKS 



In this paper we have performed a BCMA study of the magnetically catalyzed fermion and scalar masses in a 
(3+l)-dimensional Abelian Higgs- Yukawa theory in the presence of a constant magnetic field. Our results show that 
even in this multiple-coupling theory, the discrepancy between the masses obtained within the CMA and within the 
more accurate BCMA is not very significant, being the difference in the mass square of at most one order of magnitude. 
We find that the region where CMA and BCMA results exactly coincide is defined by the condition t = vln(^) <C 1. 

The BCMA calculations led to fermion masses many orders of magnitude larger than those obtained in the QED case, 
thereby confirming, within a more accurate approximation, that the Yukawa interactions strengthen the generation 
of the dynamical fermion mass by several orders of magnitude, a claim done in previous papers [12,13] based only on 
CMA results. 

As mentioned in the Introduction, a motivation for the inclusion of fcrmions-scalar interactions in the study of the 
magnetic catalysis was to find out if this phenomenon could influence the phenomenology of the early universe. A 
fundamental question here to understand is whether the strengthening of the mass by the fermion-scalar interactions 
may have any impact in the electroweak phase transition. For this effect to be of any significance for the electroweak 
physics, a condition has to be met: during the electroweak transition the universe has to be permeated by a primordial 
magnetic field strong enough as to induce, even at temperatures comparable to the electroweak critical temperature, 
a modification in the value of the fermion mass. 

We should keep in mind that at temperatures below, but close enough, the critical temperature for the electroweak 
spontaneous symmetry breaking, the fermion masses generated through the Higgs mechanism are very small, since the 
transition is expected to be either second order or weakly first order. Then, if the magnetic field is much larger than 
these tiny masses, the fermions will be mainly constrained to their LLL and the MC can be fully operative. However, 
this is only true if the thermal fluctuations are not as important as to take the fermions out of the LLL. Another way 
to put this is to say that the critical temperature at which the magnetically induced fermion mass evaporates has to 
be larger than the electroweak critical temperature. 

Magnetic fields may have well been present at the early universe. In fact, there are very plausible arguments favoring 
the existence of primordial magnetic fields that can serve as the source of the seed fields required to explain the observed 
magnetic fields in galaxies and clusters of galaxies [23] . The literature on this topic is rich in possible primordial fields 
generating mechanisms, and many of them can produce very strong fields at and before the electroweak transition 
[36,37]. 



4 The agreement between CMA and BCMA results in QED have been also proved using an improved ladder approximation 
of the gap equation [9], on which the one-loop photon propagator is used in the gap equation, instead of the bare photon 
propagator. 
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Although the model used in our calculations lacks the complexity of the electroweak theory, it shares some common 
features with the electromagnetic sector of the electroweak model, and as so we expect that any conclusion drawn 
within our model can be seen as an indication (even if qualitative) of what the relevance of the effect would be in the 
electroweak context. 

Taking into account that the critical temperature for the vanishing of the magnetically catalyzed fcrmion mass is 
typically of the order of the value of the dynamical mass at zero temperature [7,12], that is T ~ rnaiT = 0), and that 
a reasonable estimate [36,37] for the primordial magnetic field at the electroweak scale is ~ 10 24 G, one obtains, for the 
values of X y and A that gives rise to the largest zero-temperature dynamical mass, that T c ~ IGeV <C T ew ~ lOOGeV. 
Hence, no magnetically induced mass would be present at the electroweak temperature because temperature effects 
override field effects at this scale. Unless new sources of extremely large B » T 2 primordial magnetic fields can be 
identified in the future, these results indicate that the MC has no relevance during the electroweak transition. 

Nevertheless, the outcomes of this work may be important for applications of the HY model in situations where 
magnetic field effects are present at sufficiently low temperatures. We expect that they will be particularly relevant 
in condensed matter applications. As mentioned in the Introduction, a HY theory has been proposed [22,26] to 
describe the observed emergence of a secondary quasiparticle gap in high-T c superconductors at certain doping levels. 
According to recent experiments [27], the secondary gap can be also triggered by an applied magnetic field. The 
resemblance of this behavior with the MC is intriguing and deserve a thorough investigation. Such an study, in turn, 
will require the extension of the results of the present paper to the two-dimensional case in order to make quantitative 
predictions that can be compared with the experiment. 
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